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Abstract 

Based on the multiplicative rule of Schubert classes obtained in [Dua], 
we present an algorithm computing the product of two arbitrary Schubert 
classes in a flag variety G/H^ where G is a compact connected Lie group 
and H d G IS the centralizer of a one-parameter subgroup in G. 

Since all Schubert classes on G/H constitute an basis for the integral 
cohomology H*(G/H), the algorithm gives a method to compute the co- 
homology ring H*{G/H) independent of the classical spectral sequence 
method due to Leray [Li,L2] and Borel [Boi, B02]. 

2000 Mathematical Subject Classification: 14N15; 14M10 (55N33; 22E60). 
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1 Introduction 

This paper presents an algorithm computing the integral cohomology ring of a 
flag manifold G/H, where G is a compact connected Lie group and H (Z G is the 
centralizer of a one-parameter subgroup. 

The determination of the integral cohomology of a topological space is a 
classical problem in algebraic topology. However, since a flag manifold G/H 
is canonically an algebraic variety whose Chow ring is isomorphic to the integral 
cohomology H*{G/H), a complete description for the ring H*{G/H) is also of 
fundamental importance to the algebraic intersection theory of G/H ([K,S2]). 

In general, an entire account for the integral cohomology H*{X) of a space 
X leads to two inquiries. 
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Problem A. Specify an additive basis for the graded abelian group 
H*{X) that encodes the geometric formation of X (e.g. a cell decom- 
position of X) . 

Problem B. Determine the table of multiplications between these 
base elements. 

It is plausible that if X is a flag manifold G/H, a uniform solution to Problem 
A is afforded by the Basis Theorem from the Schubert enumerative calculus [S2]. 
This was originated by Ehresmann for the Grassmannians Gn^k of fc-dimensional 
subspaces in C" in 1934 [E] , extended to the case where G is a matrix group by 
Bruhat in 1954, and completed for all compact connected Lie groups by Chevalley 
in 1958 [Ch2]. We briefly recall the result. 

Let W and W' be the Weyl groups of G and H respectively. The set W/W' 
of left cosets of W' in W can be identifled with the subset of W: 

W = {w eW \ l{wi) > l{w) for all Wi e wW'}, 

where I : W ^ Tj is the length function relative to a flxed maximal torus T in G 
[BGG, 5.1. Proposition]. The key fact is that the space G/H admits a canonical 
decomposition into cells indexed by elements of W 

(1.1) G/H= U_X^, dimXy,^2l(w), 

with each cell the closure of an algebraic affine space, known as a Schubert 
variety in G/H [Chs, BGG]. 

Since only even dimensional cells are involved in the decomposition (1.1), the 
set of fundamental classes [X^j] e H2i{w){G/H), w e W, forms an additive basis 
of H,{G/H). The cocycle class e H^^^'^^G/H), w e W, defined by the 
Kronecker pairing as (P^, [X„]) = 6w,u, w,u & W, is called the Schubert class 
corresponding to w. The solution to Problem A can be stated in (cf. [BGG]) 

Bcisis Theorem. The set of Schubert classes {P^ \ w e W} consti- 
tutes an additive basis for the ring H*{G/H). 

One of the direct consequences of the basis Theorem is that the product of 
two arbitrary Schubert classes can be expressed in terms of Schubert classes. 
Precisely, given u,v E W, one has the expression 

Pu-Pv^ E _<,vP^,<,v^^ 

l{w)=l{u)+l{v),weW 

in H*{G/H). Thus, in the case of X — G/H, Problem B has a concrete form. 
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Problem B'. Determine the structure constants a^^ of the ring 
H*{G/H) for w,u,veW with l{w) = l{u) + l{v). 

Originated in tiie pioneer works of Schubert on enumerative geometry from 
1874 and spurred by Hilbert's fifteenth problem, the study of Problem B' has 
a long and outstanding history even for the very special case G = U{n) and 
H = U{k) X U{n — k), where U{n) is the unitary group of rank n (cf. [K]). The 
corresponding flag manifold G/H is the Grassmannian Gn,k of A;-planes through 
the origin in C", and the solution to Problem B' is given by the classical Fieri 
formula^ and the Littlewood-Richardson rule"^. We refer to the articles [KL] by 
Kleiman-Laksov and [St] by Stanley for full expositions of these results respec- 
tively from geometric approach and from combinatorial view-point. 

During the past half century, many achievements have been made in extending 
the knowledge on the a^^ from the Gn.k to flag manifolds of other types. See [Chi] , 
[Mo], [BGG], [D], [LS2],' [HB], [KK], [Wi], [BS1-BS3], [S2], [PRi- PR3], [Bi]. 

Early in 1953, Borel introduced a method to compute the cohomology algebra 
H*{G/ H;M.) (with real coefficients) using spectral sequence technique [Boi, B02, 
B, TW, W]. In the results so obtained the algebra H*{G/H;'R) was character- 
ized algebraically in terms of generators-relations, in which the basis theorem 
that implies the geometric structure of the space G/H was absent^. In recent 
years, in order to recover from Borel's description of the algebra H*{G/H; M) the 
polynomial representatives of Schubert classes so that explicit computation for 
the a^^ is possible, various theories of Schubert polynomials were developed for 
the cases where G is a matrix group and if C G is a maximal torus (cf. [S2], 
[LSi], [Be], [BH], [BJS], [FK], [FS], [Fu], [LPR], [Ma]). 

Combining the ideas of the Bott-Samelson resolutions of Schubert varieties 
[BS, Han] and the enumerative formula on a twisted product of 2 spheres devel- 
oped in [DU2] , the first author obtained in [Dus] a formula expressing the structure 
constant a^,, in terms of Cartan numbers of G. It was also announced in [Dua] 

^In order to find a formula for the degrees of Scliubert varieties on the Grassmanian, Schubert 
himself developed a special case of the Fieri formula [K] . 

^Classically, the Littlewood-Richardson rule describes the multiplicative rule of Schur sym- 
metric fimctions. It was first stated by Littlewood and Richardson in 1934 [LR] and completly 
proofs appeared in the 70's (see "Note and references" in [M, p. 148]). Lesieur noticed in 1947 
[L] that the multiplicative rule of Schubert classes in the Grassmanian formly coincides with 
that of Schur functions. That is, the Littlewood-Richardson rule can also be considered as the 
rule for multiplying Schubert classes in the Grassmanians. 

^In the intersection theory, the basis Theorem is important for it guarantees that the rational 
equivalence class of a subvariety in G/H can be expressed in term of the base elements and 
therefore, the intersection multiplicities of arbitrary subvarities in G/H can be computed in 
terms of the a!f ... 
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that, based on the formula, a program to compute the be compiled. This 

paper is devoted to explain the algorithm in details. 

Consequently, the algorithm gives a method to compute the integral cohomol- 
ogy ring H*{G/H) independent of the classical spectral sequence method due to 
Leray [Li,L2] and Borel [Boi, B02]. It has also served the purpose to indicate our 
further programs computing Steenrod operations on G/H [DZ], and multiplying 
Demazure basis (resp. Grothendieck basis) in the Grothendieck cohomology of 
G/H [DU4]. 

The paper is so arranged. In Section 2 we recall the formula for the a^^ 
from [Dus]. In Section 3 we resolve Problem B' into two algorithms entitled 
^^Decompositions" and "L-H coefficients". The functions of the algorithms are 
implemented respectively in section 4 and 5. 

Explicit computation in the cohomology (i.e. the Chow ring) of such classical 
spaces as flag varieties is not only required by the effective computability of 
problems from enumerative geometry [K], but also related to many problems 
from geometry and topology [IM, H, Dui]. In order to demonstrate that our 
algorithm is effective, samples of computational results from the program are 
explained and tabulated in Section 6. 

It is worth to mention that there have been excellent codes for multiplying 
Schubert classes in the Grassmannians Gn,k (cf- the programs SYMMETRICA 
at Bayreuth 



< http: / /www mathe2.uni-bayeuth.de>; 



the programs ACE at Marne La Valle 

< pittp: / /phalanstere. univ- mlv.fr >, 



and the program LITTLWOOD-RICHARDSON CALCULATOR at Aarhus 

< http://home.imf.au. dk/abu'di] >). 



Instead of being type-specific, our program applies uniformly to all G/H. 



2 The formula 

This section recalls the formula for the a^^ from [DU3]. A few preliminary nota- 
tions will be needed. Throughout this paper G is a compact connected Lie group 
with a fixed maximal torus T. We set n = dimT. 
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Equip the Lie algebra L{G) of G with an inner product (, ) so that the adjoint 
representation acts as isometrics of L{G). The Cartan subalgebra of G is the 
Euchdean subspace L(T) of L(G') [Hu, p.80]. 

The restriction of the exponential map exp : L{G) — ^ G to L(T) defines a 
set D{G) of m = |(dimG — n) hyperplanes in L{T), i.e. the set of singular 
hypcrplanes through the origin in L{T). These planes divide L{T) into finitely 
many convex cones, called the Weyl chambers of G. The refiections a of L{T) in 
the these planes generate the Weyl group W of G. 

Fix, once and for all, a regular point a e L(T)\ U L and let A = 

LeD{G) 

■ ■ ■ , f3n} be the set of simple roots relative to a [Hu, p. 47]. For a 1 < i < n, 
write ai E W for the rcfiection of L{T) in the singular plane Lf^. G D{G) corre- 
sponding to the root jSi. The (Xj are called simple reflections [Hu, 42]. 

Recall that for 1 < i, j < n, the Cartan number /?j o Pj =: 2{Pi, Pj)/{/3j, /3j) is 
always an integer (only 0, ±1, ±2, ±3 can occur) [Hu, p. 39, p. 55]. 

It is known that the set of simple refiections {cTj \ I < i < n} generates W. 
That is, any w E W admits a factorization of the form 

(2.1) W = (Jj, O • • • O (7j^, . 

Definition 1. The length l(w) of a w e is the least number of fac- 
tors in all decompositions of w in the form (2.1). The decomposition 
(2.1) is said reduced if k — l{w). 

If (2.1) is a reduced decomposition, the k x k (strictly upper triangu- 
lar) matrix A^, = {as,t) with 

_ r if s > t; 

is called the Cartan matrix of w associated to the decomposition 
(2.1). 

Let Z[a;i, • • • ,Xk] = (Br>o'^[xi, ■ ■ ■ ,Xk]^^^ be the ring of integral polynomials 
in xi, • • • , Xk, graded by | |= 1. 

Definition 2. Given an kxk strictly upper triangular integer matrix 
A = (ttij), the triangular operator associated to A is the homomor- 
phism Ta : Z[a;i, ■ ■ • ,0;^]*^'^^ — >• Z defined recursively by the following 
elimination laws. 

1) a h eZ[xi,- ■ ■,Xk-i]^''\ then Ta(/i) =0; 

2) if = 1 (consequently A = (0)), then Ta{xi) = 1; 

3) a he Z[xi, Xk-i]^''-'''^ with r > 1, then 

TA{hxl) = TA'{h{ai^kXi H h ak-i^kXk-iY~^), 
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where A' is the {{k — 1) x {k — 1) strictly upper triangular) matrix 
obtained from A by deleting the k^^ column and the k*^ row. 

By additivity, Ta is defined for every / G Z[a;i, • • • ,Xk\^''^ using the unique 
expansion / = T,hrXl with hr G Z[xi, ■ ■ ■ , Xk-i\^''~^^ ■ 

Example. Definition 2 implies an effective algorithm to evaluate Ta- 
For k — 2 and >li = ^ ^ ^ ^ , then T^^ : Z[a;i, X2]^'^^ ^ Z is given by 

Ta,{xI)=0, 

Ta^{xiX2) = Ta[{xi) = 1 and 
^A^ixl) = TA'^{axi) = a. 



a b 

For A; = 3 and ^2 = I c , then ^2 = and : a;2, Xs]*^^) 







Z is given by 



r r^n^ra rax _ r 0, if = and 

IaA^I 2^2 2^3 J - I TAAxVxl'ibx, + CX2Y'-^), if > 1, 

where ri + r2 + = 3, and where Ta^ is calculated in the above. 



Assume that w = o • • • o cTjj. , 1 < ii , • • • ,ik <n, is a reduced decomposition 
of w G W, and let = {as,t)kxk be the associated Cartan matrix. For a subset 
= bi' ■ ■ ■ ^ [1) ■ ■ ■ ) ^] we put I L 1= r and set 

ctl = o • • • o (7j^.^ eW; Xl^ Xj^ ■ ■ ■ Xj^ G Z[xi, • • • , x^]. 

The solution to Problem B' is (cf. [Dus]) 

The formula. If u,v &W with l{w) — l{u) + l{v), then 
<,v-TaM E xl){ E xk)], 

\L\=l{u), aL=u \K\=l{v), aK=v 

where L,K C [1, • • • ,k]. 



The subsequent sections are devoted to clarify the algorithm implicitly con- 
tained in the formula. 
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3 The structure of the algorithm 

Let L{T) be the Cartan subalgebra of G and let A = {f3i, ■ ■ ■ , f3n} C L{T) be 
the set of simple roots of G relative to the regular point a G L{T) (cf. Section 
2). The Cartan matrix of G is the nxn integral matrix C = {cij)nxn defined by 

It is well known that (cf. [Hu]) 

Fact 1. All simply comiccted compact semi-simple Lie groups are 
classified by their Cartan matrices. 

For a subset K = [ii, ■ ■ ■ , C [1, ■ ■ ■ , n] let 6 e L(T)\{0} be a point lying 
exactly in the singular hyperplanes L^.^ ,•''■> Lj3i^ ; namely, 

(3.1) b e n U (e L{T)\ U L^, iiK = iD) 

where J is the complement of i^' in [1, ■ ■ ■ Denote by Hk the centralizer of 
the 1-parameter subgroup {exp(t6) | t e R} in G. It can be shown that (cf. [BHi, 
13.5-13.6])) 

Fact 2. The isomorphism type of the Lie group Hx depends only on 
the subset K and not on a specific choice of b in (3.1). Further, every 
centralizer H of a one-parameter subgroup in G is conjugate in G to 
one of the subgroups Hk. 

By Fact 2 we may assume that H is of the form Hk for some K G [1, • • • ,n]. 
Summarizing Fact 1 and 2 we have 

Lemma 1. A complete set of numerical invariants required to deter- 
mine a flag manifold G/H consists of 

1) a Cartan matrix G — (cjj)„xn (to determine G); 

2) a subset K — [ii, • • • ^i^] C [1, • • • ,n] (to specify H C G). 

The implementation of our program essentially consists of two algorithms, 
whose functions may be briefed as follows. 

Algorithm A. Decompositions. 

Input: A Cartan matrix C = {cij)nxn, and a subset K G [1, . . . ,n\. 
Output: The coset W being presented by a reduced decomposition for each 

w eW. 



7 



Remcirk 1. In [Ste, Section 1] Stembridge described an algorithm for the 
problem of finding a reduced decomposition for a given w & W. This requests 
less than what Algorithm A concerns. 

Algorithm B. L-R coefficients. 

Input: u, V, w E W witti l{u) + l{v) — l{w). 

Output: a"^^^ e Z. 

The details of these algorithms will be given respectively in the coming two 
sections. 

It is clear from the above discussion that our algorithms reduce the structure 
constants a^.^ directly to the Cartan matrix C = (Qj)„xn and the subset K C 
[1, . . . , n]: the simplest and minimum set of constants by wfiidi all flag manifolds 
G/H are classified (cf. Lemma 1). Because of this feature it is functional equally 
for computations in all G/H. 



4 Algorithm A 

We show in 4.1 the fashion by which the Weyl groups W CW arise from the 
Cartan matrix C = {cij)nxn and the subset K G [I, ■ ■ ■ ,n]. In 4.2 a numerical 
representation for W is introduced. Based on the terminologies developed in 4.1 
and 4.2, Algorithm A is given in 4.3. 

4.1. Constructing the Weyl groups W' C W from the Cartan matrix. 

Let r be the free Z- module with n generators cui, - ■ ■ ,u!n, and let Aut{r) be the 
group of automorphisms of F. 

Given a Cartan matrix C = {cij)nxn of a Lie group G with rank n, define n 
endomorphisms of F (in term of Cartan numbers) by 

(4.1) a,{u,) = I if M ^; . , 1 < ^ < n. 

It is straightforward to verify that cr^ = Id. In particular, CTfe G Aut{T). 

Lemma 2. The subgroup of Aut{r) generated by ai,- ■ • , an is iso- 
morphic to W, the Weyl group of G. 

For a subset K C [1, ■ ■ ■ the subgroup W' of W generated by 
{(Ti I i e K} is isomorphic the Weyl group of Hk (cf. Section 3). 

Proof (cf. Proof of Theorem 1 in [DZZ]). Let t be the real vector space 
spanned by a;i, • • • ,ujn\ namely, i = F (8) IR. In term of the Cartan matrix C — 
{cij)nxn we iutroduce in t the vectors • • • , /3„ by 
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A — Cil^l + • • • + CinUJn, 



and define an EucHdean metric on t by 

(4.2) 2(A,(^) = Q,; (A,/3i) = l. 

Then 

(a) t can be identified with the Cartan subalgebra L{T) of G under which the 
vectors • • • , corresponds to the set A of simple roots of G (cf. Section 2); 

(b) with respect to the metric (4.2), the induced action of ak on t — L{T) is 
the reflection in the hyperplane L^^ perpendicular to the f3k', 

(c) under the identification t = L(T) specified in (a), the basis Ui, ■ ■ ■ ,cj„ of 
r agrees with the set of the fundamental dominant weights relative to A [Hu, 
p. 67] (Geometrically, positive multiples of wi, • • • ,a;„ form the edges of the Weyl 
chamber in L{T) corresponding to A). 

Lemma 2 follows directly from (b) and (c).n 

4.2. A numerical representation of Weyl groups In the theory of Lie 
algebras the vector S — u;i-\ ha^neTct^r^lRis well known as a strongly- 
dominant weight [Hu, p. 30]. For aw &W consider the expression in F 

w{5) = biUi H h bntOn: h G Z. 

Definition 3. The correspondence 6 : ^ by h{w) — {h-i, ■ ■ ■ ,hn) 
will be called the numerical representation of W. 

Lemma 3. The numerical representation h -.W ^ is faithful and 
satisfies bi ^ for all w & W and 1 < i < n. 

Proof. By (c) in the proof of Lemma 2, 5 G t is a regular point in the 

Weyl chamber determined by A. Lemma 3 comes from the geometric fact that 
the action of the Weyl group W on the orbit of any regular point is simply 
transitive. □ 

The formula (4.1), together with additivity of the ak, is sufficient to compute 
the coordinates of b{w) from the Cartan numbers and any decomposition of G 
W into products of the cxj, as the following algorithm shows. 

Algorithm 1. Computing b{w). 

Input: A sequence 1 < ii, ■ ■ ■ ,im ^ n-. 

Output: b{w) for w = (Tij o • • • o (ji^. 

Procedure: Begin with the sum pq = ui + • • • + cUn- 

Step 1. Substituting in po the term by cui^ - T,i<j<nCi^j(^j to get pi, 

Step 2. Substituting in pi the term a;^^,^ by u}i^_i - '^i<j<nCim-ij<^j to get p2; 
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Step m. Substituting in Pm-i the term cui^ by cUjj — Tji<j<nCiijUJj to get Pm] 
Step m+1. If Pm — bicui + • • • + bnCUn then b{w) — {hi, ■ ■ • , 6„). 

Wc conchidc this subsection with two useful properties of the numerical rep- 
resentation of a Weyl group given in Definition 3. Let I : W ^ Z he the length 
function on 1^. As in Section 1 we identify W with the subset of W 

W^{weW\ l{w) < l{u) for all u e wW'}. 

Lemma 4. Let w &W he with h{w) = [hi, ■ ■■ ,hn) and b{w~^) ~ 
{h,--- X)- Then 

(i) l{aiw) — l{w) — 1 if and only if hi < 0; 
(a) w & W if and only if bi > for all i e K. 

Proof. The metric on L{T) yields the relations 

(4.3) (a;,,/?,V(/?,-,/?,)) = (5,,- 

between the simple roots j3j and the corresponding fundamental dominant weights 
cui [Hu,p.67]. By [BGG, 2.3 Corollary], /((Tjw) = if and only if A) < 

0. The latter is equivalent to 6j < in views of (4.3) and w{5) = hiuj\ + - ■ ■ + hnUJn- 
This verifies (i). 

Similarly, assertion (ii) follows from the following alternative description for 
W (cf. [EGG, 5.1. Proposition, (iii)]) 

W^{weW\ (w-i((5), A) > for aiJ i e K}n 



4.3. Construction of the coset W = W/W. Let / : ^ Z be the length 
function on W. We put w'' = {w e W \ l{w) = /c}, /c = 0, 1, 2, ■ ■ ■ . Then, as is 
clear, = ]J . The problem concerned by Algorithm A may be reduced to 

fe>0 

Problem C. Enumerate elements in (i.e. in W), k > Q, by their 
reduced decompositions. 

Before presenting Algorithm A (i.e. the solution to Problem C) we note that 

(4.4) If the set w'^ is given in term of certain reduced decompositions of its 
elements, then W becomes an ordered set with the order specified by 

if there exists some s < k such that it = jt for alH < s but ig < js- 
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(4.5) If X and Y are two ordered sets, then the product X xY is furnished with 
the canonical order as: 

"(x, y) < {x', y') if and only ii x < x' or x — x' but y < y'" . 

The solution for Problem C is known when k = 0,1 

= {id}; W' = {aj I j e J}, 

where id is the identity of W and where J is the complement of in [!,••• ,n]. 

fc— 1 

In general, Algorithm A enables one to build up W from W 

Algorithm A. Decompositions. 

fc— 1 

Input. The set W being presented by certain reduced decompositions of 
its elements. 

Output. The set W being presented by certain reduced decompositions of 
its elements. 

fc— 1 

Procedure: Set V — {1, ■ • ■ ,n} x W . Repeat the following steps for all 
elements in V in accordance with the order on V (cf. (4.5)). Begin with empty 

sets 5 = 0, i? = 0. 

Step 1. For a. v = (i, cTj^ o • • • o CTj^ J e V form the product w = (Xj o cjjj o • • • o 
Step 2. Call Algorithm 1 to obtain 

b{w) = {bi, ■■■ ,bn) and b{w~^) = (6i, • • • ,6„); 

Step 3. If 1) bi_< 0; 

2) 6i > for all i E K; 

3) (6i,---,6„)^i?, 

add (jj o CTj^ o ■ ■ ■ o o'jj._^ to S; add b{w) = (6i, • • • , 6„) to R; 



The program terminates at S = W . 

Explanation. We verify the last clause in Algorithm A. Firstly, Lemma 7 in 

[DZZ] claims that any w E W admits a decomposition w = ai o ai^^ o ■ ■ ■ o ai^_^ 
for some (i, cij^ o ■ ■ ■ o (Tj^. J G V . This explains the role the set V plays in the 
algorithm. Next, the first two conditions in Step 3 guarantees that a^o Ui-^o ■ ■ ■ o 

G W by Lemma 4. Finally, the third constraint in Step 3 rejects a second 

reduced decomposition of some w G being included in (by Lemma 3).n 

Remark 2. li K = 0, then W = W (the whole group). In this case H = T 
(a maximal torus in G) and Step 2 and 3 in Algorithm A can be simplified as 
Step 2. Call Algorithm 1 to obtain b{w) = (6i, • • • ,bn); 
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Step 3. If 6j < and if (&i, •••,&«) ^ R, add aj o o • • • o aii^_^ to S; add 
= (61, • • • to R. 

Remark 3. Based on the word representation of Weyl gronps, a different 
program solving Problem C was given in [DZZ]. In comparison, the use of the 
numerical representation simplifies the presentation of Algorithm A. 

5 Algorithm B 

Algorithm A presents us the coset = U by certain reduced decomposi- 

A;>0 

tion of its elements. Based on this we explain L-R coefficients, the algorithm 
computing a^„. 

By the notion L C [!,••• , A;] , |L| = r, we mean that L is a sequence (ji, • • • , jV) 
of r integers satisfying 

1 < ii < • • • < ir < A;. 

For two integers 1 < r < A; let the set V{k, r) = {L | L C [1, • • • , /c], = r} be 
equipped with the obvious ordering (cf. (4.3)). 

For di, w — Ui^o ■ ■ ■ o Gii^ &W and a e 1^ with l{u) = r < k, we set 

(5.1) Pw{u)^ E XL eZ[xi,--- ,Xk]^''\ 

L£V{k,r),aL=u 

where o"l = <^ij-^ o ■ ■ ■ o if L = [71, • • " tJA- Using these notations our formula 
(cf. Section 2) can be simplified as 

(5.2) al,^TAMu)pM]. 

We begin by pointing out that (5.1) suggests the following algorithm specify- 
ing the polynomial Pw{u). 

Algorithm 2. Computing Pw{u) E Z[a;i, ■ ■ ■ ,Xk]^^\ 

Input: -u; = (Tjj o • • • o cTj^ E W and u G with b{u) = {bi, ■ ■ ■ ,hn). 

Output: Pw{u). 

Procedure: Repeat the following steps for all L e V{k,r) in accordance with 
the order on V{k, r). Initiate the polynomial p — p{xi, ■ ■ ■ , Xk) as zero. 

Step 1. For a L G V{k,r) call algorithm 1 to get b{aL)] 
Step 2. If b^ai) = b(u) add Xl to p. 
The program terminates at p = Pw{u). 
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If A — {aij)kxk is matrix of rank k and ifl<r<A; — 1, then the notion 
{aij)rxr clearly stands for the matrix of rank r obtained from A by deleting the 
last {k — r) rows and columns. 

Let A = {aij)kxk be a strictly upper triangular integral matrix of rank k. 
Consider the triangular operator Ta : Z[xi, ■ ■ ■ , Xk]^''^ —>■ 'Z^ given in Definition 2. 

Algorithm 3. Computing Ta : Z,[a;i, • • • ,0;^]^'^) Z. 

Input: A strictly upper triangular integral matrix A — {aij)kxk and a poly- 
nomial p = p{xi, ■ ■ ■ , Xk) e • • • , Xk]^^"^ 
Output: Ta{p) e Z. 
Procedure: Recursion. 

Step 1. Express p as a polynomial in Xk] i.e. 

P = ho + hiXk + hrXl, hr e Z[Xi, ■■■ , Xk-lY^~''\ 

2<r<k 

and set 

Pi = ^1 + E K{ai,kXi H h ak-i,kXk-iY~^{& • • • , Xk-i]'^''''^^). 

2<r<k 

Step 2. Repeat step 1 for Ai — (aij)(fc-i)x(fe-i) and p = pi to get p2 G 

Step k+1. li pk — a & Z, then T^(p) = a. 

Algorithm B. L-R coefRcients. 

Input: w — Gi^o ■ ■ ■ o e , {u, v) e x l^'^ ^ 

Output: a^^^, G Z. 

Procedure: Let be the Cartan matrix of w related to the decomposition 
(it can be read directly from the Cartan matrix of G and the decomposition 
w — Gi^o ■ ■ ■ o Gi^. cf. Definition 1). 

Step 1. Call algorithm 2 to get Pw{u) and 

Step 2. Call algorithm 3 to get Ta^{Pw{u) -Pwiv)). 

Step 3. If Ta^{Pw{u) -Pwiv)) = a, then a;^„ = a (by (5.2)). 

Remcirk 4. Based on Algorithm B, a parallel program to expand the product 

P ■ P = o'^ P 

for given {u, v) e x w'' ^ can be easily implemented. The order on w'' can 
be employed to assign each w eW a computing unit. 
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6 Computational examples 



Our algorithm is ready to apply to computation in flag manifolds. Recall that 
all compact connected semi-simple irreducible Lie groups fall into four infinite 
sequences of matrix groups 

SU{n)] S0{2n)] 50(2n+l); Sp{n), 

as well as the five exceptional ones 

G2, -F4, Eq, Et, Es- 

The flag manifolds associated to matrix groups have been studied extensively 
during the past decades. Here, we choose to work with certain flag manifolds 
related to the exceptional Lie groups E^, n — 6,7, 8. 

Fix a maximal torus T" C En and let W{n) be the Wcyl group of En. Then 

2^3^5 if n = 6; 78 if n = 6; 

(6.1) \W{n)\ = { if n = 7; dim^ E^ = { 133 if n = 7; 

2^^35527 if n = 8; 248 if n = 8, 

where \A\ stands for the cardinality of the set A. Assume that the set of simple 
roots A = • • ■ , Pn} of En is given and ordered as the vertices of the Dynkin 
diagram of En pictured in [Hu, p. 58], and let K C {1,2, •• • ,n} be the subset 
whose complement is {2}. We have the following information on the subgroup 
group Hk- 

(a) the semisimple part of the subgroup Hk C En is SU{n), the 
special unitary group of order n; 

(b) Hk admits a factorization into the semi-product Hk — S^-SU{ri), 
where is a circle subgroup of the maximal torus in En; 

(c) if W'{n) C W{n) is the Weyl group of Hk, then \W'{n)\ = n\. 
Consequently, if one write W{n) for the coset W' (n) in W{n), one has 

_ 2^32 if n = 6; 42 if n = 6; 

(6.2) \W{n)\ = { 2632 if n = 7; dimR En/HK = { 84 if n = 7; 

2^335 if n = 8. 194 if n = 8. 

Geometrically, W{n) parameterizes Schubert classes on En/HK (i.e. the Basis 
Theorem) . 

The subset of W{n) consisting of elements with length r is denoted W\n). By 
(4.4), if the W (n) is presented by its elements each with a reduced decomposi- 
tions, then it naturally becomes an ordered set and therefore, can be alternatively 
presented as 
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(6.3) W\n) = {wr,i\l <i <\W"\}. 

In table A„ below, we present elements of W{n) with length r < 10 both in 
terms of their reduced decompositions produced by Algorithm A, and the index 
system (6.3) imposed by the decompositions. 

The index (6.3) on W (n) is useful in simplifying the presentation of the 
intersection multiplicities a^„. By resorting to this index system we list in table 
B„ ( n = 6, 7, 8) all the a^,, with l{w) = 9 and 10 produced by Algorithm B. 

Table Ag. Reduced decomposition of elements in W{6) with length< 10 





decomposition 






decomposition 








^7,5 


£T5£T4a3£r6£75£74£7'2 


W2,l 


Cr4(J2 




^8,1 


£7lCr2£74£73£76a5£74£r2 




(73(74(72 




^^8,2 


£7l £75£74£73f76£75£74(72 


W3,2 






W^8,3 


0"2 C4C3C5C4C2 


W4,l 


C1C3C4C2 






£T2£75£74£73£76£75£74£72 


U!4,2 


(7,3(75(74(72 




Ws,5 


£T3£ri£74£T3£Tg£T5£T4£T2 


W4,3 


(76C75C74a2 




W9,l 


£7l£72£75£74£T3a6£T5£74£72 


^5,1 


(71(73(75(74(72 




W9,2 


£72 C3 C 1 (74 (73 £76 (75 (74 (72 


^5,2 


C3C6C5C4C2 




W9,3 


£73 £7i£75£74(73(76(75(74(72 




£74 0-3 0-5 £74 £72 




W9,4 


£T4£72£73£7i£74£73£75£74£72 




£Tl£T3£T60"5£74£72 




^9,5 


£74 £72 £75 £74 £73 £70 £75 £74 £72 


W6,2 


£Ti£T4£T3£75£74£T2 




^10,1 


£7 1 £74 £72 £75 £74 £73 £76 £75 £74 £72 


^6,3 


£72(74(73£75£74(72 




^10,2 


£72 £73 £7 1 £75 (74 £73 £76 (75 (74 £72 


W6,4 


C4C3f76£75£74£72 




^10,3 


O'S C4C2 C5C4(73(76(75£74£72 


^7,1 


£Tl£T2C74£T3£T5£T4£T2 




^•'''10,4 


£74£72£73£7i£74£73£76£75£74£72 


W7,2 


£7lC74C73£76a5£T4£r2 




W^10,5 


£74£73£7i£75£74£73£76£75£74£72 


W7,3 


£12 £7,1 £73 £75 £75 CT4 £72 




^•'''10,6 


£75 £74 £72 £73 £7^ £74 £73 £75 £74 £72 


U)74 


aj,a 1(7. i_o':iO'[,o' 1(72 
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Table A7. Reduced decomposition of elements in W{7) with length< 10 







WiJ 


ClcCUllipUolllUli 






(_lfc;CUllip(Jk5lllUli 


Wl^l 


C2 


'^^7,3 


(Tl(T4(T3(T6Cr5Cr4(T2 




^9,6 


•73 (7 1 (74 £73 (77 (76 £75 (74 £72 


W2,l 


C4Cr2 


^7,4 


(T2Cr4(T3(Tg(T5(T4(T2 






0'3(7l£75£74£73£76£75£74£72 


^3,1 


CT3(T4C72 




0'3ClC4C3C5<74C2 




''^9,8 


£74 £72 £73 £7 1 £74 £73 £75 £74 £72 


^3,2 


(T5Cr4(T2 


^7,6 


Cr4(73(77(T6Cr5(T4(T2 




^^9,9 


£74 (72 (75 £74 (73 £76 £75 (74 £72 


W4^l 


(Tl(73(74(72 


W-jj 


(T5(74(73(T60'5(T4(T2 




^9,10 


(76 (75 £74 £73 (77 £76 £75 (74 £72 


W4^2 


(T3CT5(T4a2 




(7l(72(T4C73(T5(J5(T4(J2 




^10,1 


(7 1 £72 £75 £74 £73 £77 £76 (75 (74 £72 


W4,3 


C6<75C74C72 


^8,2 


C 1 (74 (T3 £77 (Tg CJ5 (74 (7 2 




U^10,2 


£7 1 £74 £72 £75 £74 £73 £76 £75 £74 £72 


W5^1 


cri(73(T5Cr4(72 


^8,3 


(T 1 (T5 (74 (T3 (Tg (75 (T4 (72 




^10,3 


•7 1 (76 £75 (74 (73 £77 £76 (75 £74 £72 


1^5,2 


^3C6C5f''4C2 


^^8,4 


'72'73C7i (74(73(75(740"2 




^^10, 4 


(72(73£7l £74£73£77£76£75£74£72 


■"^5,3 


(T4Cr3C75(T4C72 


''^8,5 


C72 C74 (73 (77 (76 (75 £74 (72 






(72£73£7i £75£74£73£76£75£74£72 


■"^5,4 


(T7(T6(75Cr4(72 


^8,6 


(72 <75 (74 (73 (76 (75 (74 (72 




■"^10,6 


£72 (76 (75 £74 (73 (77 £76 (75 £74 £72 


W6,l 


ClCr3Cr60'5C^4C2 




(73(71 (74(73(76(75<74(72 




mo,7 


£73 (71 (75 (74 (73 (77 £76 (75 £74 £72 


Wq^2 


(Ji(T4(T3Cr5Cr4(T2 


W^8,8 


C75£74(73£77£76(75£74(72 






£73£74£72£75£74£73£76(75(74£72 


^i^6,3 


<72<740'3<750'4i72 


W9,l 


(7i(72(74C73C77(76(75C74(72 






(74(72£73£7i£74£73£76£75£74£72 




C30'70'6C50'4C2 


W9,2 


(71(72(75(74(73(76(75(74(72 




W^ICIO 


£74(72£75(74(73£77£76(75£74£72 


we,5 


cr4Cr3Cr6(J5Cr4Cr2 


W9,3 


£7 1 £75 (74 (73 £77 (76 £75 (74 (72 




^^10,11 


£74£73£7i £75£74£73£76£75£74£72 




Cri(J2(T4C3<75C740"2 


Wg^4: 


£72£73(7i£74£73(76£75(74(72 




^t^l0,12 


(75 (74 £72 £73 £7 1 £74 £73 £75 £74 £72 


W7,2 


Cri(T3(T7Cr6(75(74(T2 


W9,5 


£72(75(74(73(77(76(75(74(72 









Table A.^. Reduced decomposition of elements in W{8) with length< 10. 





decomposition 




WiJ 


decomposition 






decomposition 


Wl^l 


0-2 




Wjfi 


£73£78£77£76£75£74£72 




W9,10 


£T4£72£73£7i £74a3£T5£74a2 


W2,l 


(74(72 




W7,7 


£74£73£77£76£75£74(72 




W9,ll 


(74(72(75(74(73(76(75(74(72 


W3,l 


(73£74£72 




W7,S 


£75£74£73£76£75£74£72 




W9^12 


£75 £74 (73 £78 £77 (76 £75 £74 (72 


W3,2 


(75(74(72 




W8,l 


£7 1 £72 £74 £73 £76 £75 £74 £72 




W9,13 


(76 £75 (74 £73 £77 £76 £75 £74 £72 


Wi^l 


(71(73(74(72 




W8,2 


(71(73(78(77(76(75(74(72 




Wio^l 


(7 1 £72 £74 £73 £78 £77 £76 (75 (74 (72 


^4,2 


f73£75£74£72 




W8,3 


£7i£74£73£77£76£75£74(72 




Wl0^2 


£7i £72(75£74£73(77£76£75(74(72 


^^4,3 


(76(75£74£72 




W8,4 


£7i(75(74£73£76(75(74(72 




Wl0^3 


£7 1 £74 £72 £75 £74 (73 £76 (7 5 (74 (72 


WFjS 


(7i£73f75£74£72 




W8,5 


(72 £73 £7 1 £74 £73 £75 £74 £72 




Wlo^A 


(7i £75£74£73£78a7£T6a5£T4£72 


W5,2 


Cr3(76<75(74£72 




W8,6 


(72(74£73£77£76£75(74(72 




f«10,5 


Cl<76(75£74£73(77(76(75(74(72 


W5,3 


C"4(73(75£74£72 




W8,7 


£72£75£74£73£76£75£74(72 




WlOfi 


(72(73(71 £74£73£77£76£75(74(72 


W5,4: 


Cr7(76(75(74C72 




W8,8 


£73£7i C74£73C76£75£74f72 




Wio,7 


£72 £73 £7 1 £75 £74 £73 £76 £75 £74 £72 


W(i,l 


(71(73(76(75(74(72 




W8,9 


£74£73£78£77£76£75(74(72 




Ww,8 


(72(75(74(73(78(77(76(75(74(72 


W6,2 


(7l(74£73£75£74£72 




Wg^lO 


£75£74£73£77£76£75£74(72 




Wl0,9 


(72(76£75£74£73(77£76£75(74(72 


W6,3 


(72£74£73f75£74£72 




wg^l 


£7 1 (72 £74 £73 £77 £76 (7 5 (74 £72 




''i'10,10 


£73 £7 1 (74 £73 £78 (7 7 £76 (7 5 £74 £72 


W(i,4: 


(73£77£76£75£74£72 




W9,2 


£7i£72£75£74£73a6£75a4a2 




Wlo^U 


(73 (7 1 £75 £74 £73 £7 7 £76 £7 5 £74 £72 


W6,5 


Cr4(73<76(75£74(72 




W9,3 


£7 1 £74 £73 £78 £77 £76 (75 (74 £72 




■i«10,12 


<73(74(72£75(74(73(76(75(74(72 


WQfi 


(78£77£76£75£74£72 




W9,4 


(7 1 £75 £74 £73 £77 £76 £75 £74 £72 




^^10,13 


(74(72(73(71 £74£73£76£75(74(72 


W7,l 


(71(72(74(73(75(74(72 




W9,5 


£72£73£7i £74£73£76£75f74£72 




W^10,14 


(74(72(75(74(73(77(76(75(74(72 


W7,2 


(7i£73(77(76£75(74(72 




W9,6 


C2 (74 (73 £78 £77 £76 (75 (74 £72 




Ww,l5 


<74 (73 (7i (75 (74 (73 (76 (75 (74 (72 


W7,3 


£7l£74£73£76£75£74£72 




W9J 


£72£75£74£73£77£76£75(74£72 




WlO^W 


(75 £74 £72 £73 £7 1 (74 £73 £75 (74 (72 


Wj^i 


(72(74£73£76£75£74£72 




W9,8 


£73£7i £74(73£77£76£75(74£72 






£70 £75 £74 £73 £78 £77 £76 £75 £74 (72 


W7,5 


£73£7i£74£73£75£74£72 




W9,9 


£73£7i £75£74£73£76£75£74£72 
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Table Bg. L-R coefficients for Eq/S^ ■ SU{Q) 



u 


f 




w 


G iy«(6) 






Wg,2 


W9,3 








W8,l 


1 


1 













W8,2 


1 





1 










W8,3 





1 





1 







W8,4 


1 











1 




U'8,5 





1 


1 








W2,l 




1 


2 





1 





W2,l 


W7,2 


2 


2 


2 








W2,l 


U'7,3 


2 


1 








1 


W2,l 


W7A 





2 


1 


1 





W2a 


W7,5 


2 





1 





1 


W3,l 


W6,l 


1 


1 


1 










We.2 


1 


3 


2 


1 







W6,3 


2 


1 





1 







We, 4 


3 


2 


1 





1 


W3,2 


W6,l 


1 


1 


1 










we,2 


2 


3 


1 


1 







we,3 


1 


2 








1 


M'3,2 


W6,4 


3 


1 


2 





1 




W5,l 





1 


1 








Wi,l 


Wz,2 


i 


i 








u 




W5,3 


1 


1 


1 


1 





Ul4,2 


W5,i 


2 


3 


2 


1 





«^4,2 


W5,2 


3 


2 


2 





1 


Wi,2 


W5,3 


5 


5 


2 


1 


1 


Wi,Z 


W5,l 


1 


1 











W4,3 


W5,2 


1 





1 








W4,3 


W5,3 


1 


1 


1 





1 
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u 


V 




wio,i 


Wl0,2 


Wl0,3 


Wl0,4 


Wio,5 


wio.e 


Wl,l 


WgA 


1 


1 














Wl,l 


W9,2 





1 





1 








Wl,! 


Wg^3 





1 








1 





Wl,! 


W9A 











1 





1 


Wis 


W9.5 


1 





1 











W2,l 


Ws,l 


1 


2 





1 








W2,l 


U'8,2 


1 


2 








1 





W2,l 


WS,3 





1 





2 





1 


W2,l 


WSA 


2 


1 


1 











W2,l 







2 





1 


1 





W3,l 


W7A 





2 





1 





1 


^3,1 


Wt.2 


1 


3 





1 


1 





W3,l 


W7.3 


2 


1 





1 








W3,l 


W7A 





1 





2 


1 





W3,l 


^7,5 


1 


2 


1 











W3,2 


W7A 


1 


1 





2 








W3,2 


W7.2 


1 


3 





1 


1 





W3,2 


W7.3 


1 


2 


1 











W3,2 


W7A 





2 





1 





1 


W3,2 




2 


1 








1 





^4,1 


Wf,.i 





1 














W4,l 


We, 2 





1 





1 


1 





W4,l 


Wq.3 





1 











1 


W4,l 


Wqa 


1 


1 





1 








W4,2 


Wqa 


1 


2 





1 


1 





Wi,2 


We, 2 


1 


5 





3 


1 


1 


W4,2 


We, 3 


2 


2 





2 








W4,2 


We, 4 


3 


5 


1 


1 


1 





^4,3 


we,i 





1 














^4,3 


We, 2 


1 


1 





1 








^4,3 


We, 3 





1 


1 











W4,3 


We, 4 


1 


1 








1 





W5,l 


W5,l 





2 





1 


1 





W5,l 


W5,2 


1 


2 





1 








W5,l 


W5,3 


1 


3 





2 


1 


1 


W5,2 


W5,2 


1 


2 








1 





W5,2 


W5,3 


2 


3 


1 


1 


1 





105,3 


W5,3 


3 


6 





3 
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Table B7. L-R coefficients for Er/S^ ■ SU{7) 



7/ 


V 








"^9,3 


079^4 


•^9,5 


•^9,6 


"^9,7 


"^9,8 


"^9,9 


"-'9,10 


"^1,1 


n't!-. 1 

"^8,1 


1 
1 


1 
1 


n 

u 


1 


n 

u 


n 

u 


n 

u 


n 

u 


n 

u 


n 
u 


"3/ Ti 1 


<^8,2 


1 
1 


n 

u 


1 

1 


n 

u 


n 

u 


1 

1 


n 

u 


n 

u 


n 

u 


n 

u 




"^8,3 


n 
u 


1 


1 

1 


n 


n 


n 


1 

1 


n 


n 


n 




7/?o ^ 

"^8,4 


n 
u 


n 


n 


1 


n 


n 


n 


1 

1 


n 


n 


^1,1 


7/fo r 

"^8,5 


1 


n 


n 


n 


1 

1 


n 


n 


n 


n 


n 


"3/ I-i 1 
"^1,1 


"^8,6 


n 

u 


1 

1 


n 

u 


n 

u 


1 

1 


n 

u 


n 

u 


n 

u 


1 

1 


n 

u 




^8,7 


n 

u 


n 

u 


n 

u 


1 

1 


n 

u 


1 

1 


1 

1 


n 

u 


n 

u 


n 

u 




f^8,8 


n 

u 


n 
u 


1 

1 


n 


1 

1 


n 


n 


n 


n 


1 

1 


^2,1 


7/7^ 1 

"^7,1 


1 

± 


1 


n 


9 
z 


n 


n 


n 


1 

1 


n 


n 


^2,1 


0^7,2 


1 


n 
u 


1 

1 


n 


n 


1 

1 


n 

u 


n 


n 

u 


n 




'^^7,3 


9 


9 
z 


9 


9 
z 


n 


9 


9 


n 


n 
u 


n 


"^2,1 


^7,4 


9 


9 
z 


n 


1 


9 


n 


n 


n 

u 


1 
1 


n 


7/?o 1 


7/Jr7 r 


n 


n 
u 


n 


9 
z 


n 


1 

1 


1 

1 


1 

1 


n 


n 


9/ Jo 1 


^7,6 


9 


n 
u 


9 


n 


9 


1 

1 


n 


n 


n 


1 

1 


"^2,1 




n 
u 


9 
z 


9 


n 

u 


9 


n 


1 

± 


n 


1 

1 


1 

1 


^3,1 


"^6,1 


1 


1 


1 

1 


1 


n 


9 


1 

1 


n 


n 


n 


"^3,1 


^6,2 


1 

1 


1 

1 


1 

1 


Q 


n 

u 


1 
1 


9 


1 

1 


n 

u 


n 

u 


"^3,1 


<^6,3 


1 
1 


9 
z 


n 

u 


1 
1 


1 

1 


n 

u 


n 

u 


1 

1 


n 

u 


n 

u 




^^6,4 


9 
z 


n 
u 


1 


n 
u 


1 

1 


1 

1 


n 


n 


n 


n 
u 


7/Jo 1 

f^3,l 




Q 



Q 


Q 


9 
z 


9 


1 

1 


1 

1 


n 


1 

± 


1 


"^3,2 


"^6,1 


9 
z 


1 


9 


1 


n 


1 

1 


1 

1 


n 


n 


n 


7/Jo r\ 


"^6,2 


9 
z 


9 
z 


1 
1 


Q 



n 


9 


1 

1 


1 

1 


n 


n 


^3,2 


"^6,3 


9 
z 


1 
1 


n 

u 


9 
z 


1 

1 


n 
u 


n 

u 


n 

u 


1 

1 


n 

u 


t^3,2 


"-'6,4 




n 
u 
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The computations were carried out by using Mathematicae on a PC. PIV667. 
Ram 128. Win98. In general, the running time of the program depends on 



(1) the order of the cosct W; 

(2) the number of non-zero entries in the Cartan matrix of G. 

More precisely, to obtain the results in Table A„, the times consumed (in seconds) 
are 



n 


6 


7 


8 


time 


1 


1 


2 



The running times for computing all aj^^ with l{w) < 10 are respectively 
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n 


6 


7 


8 


time 


38 


115 


159 
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